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Symmetry arguments are used to describe all possible two-dimensional periodic corrugations of 
graphene ("nanomeshes") capable of inducing tangible semiconducting gap. Such nanomeshes or 
superlattices break the initial graphene translational symmetry in a way that produces mixing and 
subsequent splitting of the Dirac K and K' states. All of them have hexagonal Bravais lattice 
and are described by space groups that are subgroups of the graphene group. The first-principles 
calculations show that the gaps of about 0.5 eV can be induced at strains safely smaller than the 
graphene failure strain. 
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Graphene, a truly two-dimensional crystal of carbon 
atoms with a honeycomb lattice, exhibits giant carrier 
mobilities up to 6 x 10^ cm^/Vs [IHS], it has a record 
mechanical strength with critical deformations exceed- 
ing 20% and huge thermal conductivity reaching 1800 
WK~^m~^|4:. As such, it is very promising as a circuit 
material in electronic devices. However, the absence of a 
band gap in graphene prevents its immediate use as e.g. a 
channel in switchable devices like field effect transistors. 

Perhaps, the most straightforward way to induce a 
gap in graphene is through quantum confinement - one 
should simply cut a graphene sheet into small pieces like 
nanoribbons [3 |6] or dots [7 . However, it is difficult to 
control the band gap produced in such a way because it 
strongly depends on the size and edge geometry of the 
nanostructures. That is why many papers have been fo- 
cused on the opening up a spectral gap in a 'bulk' mate- 
rial using different kinds of periodic lattice perturbations 
with various substrates [SI |9], antidot lattices [lOVllSj . 
regular lines or areas of adsorbed foreign atoms [T1HT7] . 
external periodic electric and magnetic fields [THj, laser 
field [11], Kekule distortions [2D], etc. 

Mechanical out-of-plane sine-wave deformations cer- 
tainly belong to this class of perturbations. Recently, we 
showed that the one-dimensional (ID) periodic ripples 
can open a gap at the Dirac points only due to (i) break- 
ing of the inversion symmetry or equivalency between A 
and B sublattices and/or (ii) merging of two inequiva- 
lent Dirac points, D and —D [21\ . Breaking of inversion 
symmetry was found to have only a relatively modest 
contribution to the gap, especially when the amplitude 
A of a sine wave deformation is relatively small in com- 
parison with the period A. A tangible gap can mainly 
result from mutual annihilation of the Dirac points D 
and —D, yet this requires severe corrugations. 

Importantly, in the case of ID periodic corrugations 
the energy gap at the Fermi level cannot be induced due 
to mixing of electronic states belonging to two different 
valleys (or to two inequivalent Dirac points K and K' 
in the unperturbed graphene) even if the imposed period 
provides a momentum transfer g = K — K' coupling the 



Dirac points. The reason is that the Dirac points gener- 
ally shift from their initial positions as the amplitude of 
corrugations A increases, so that the relation g = K K' 
does not hold any longer at any finite amplitude A. 

In the present study, we show that in contrast to ID 
the 2D out of plane modulations can easily induce the 
gap at the Dirac points via K — K' mixing, provided 
that the modulations have hexagonal translational sym- 
metry. It might appear that this idea falls under the 
concept developed in Ref. ^5] of creating the gaps in 
graphene by putting the latter on top of a corrugated 
surface with a triangular landscape. In reality, however, 
there is nothing in common between the two since we 
consider the semiconducting gaps forming at the Dirac 
points due to inter-valley mixing, while the authors |22| 
studied the multiple minigaps lying above (and below) 
the Dirac energies. 

The opportunity to open up the gap through 2D 
corrugations has become especially interesting in light 
of the recent experiments showing that such corruga- 
tions ("nanomeshes") can be formed during the epi- 
taxial graphene growth on some special substrates like 
SiC(OOOl) [23H26] or transition metal surfaces [27H29] . 
When grown, for example, on Ru(OOOl), graphene forms 
a nanomesh where 25 x 25 carbon hexagons are commen- 
surate with 23 X 23 unit cells of Ru [28] . 

In what follows, we consider a graphene sheet with 
2D corrugated profile Uz{r) that forms either due to a 
lattice misfit with a substrate or follows its surface cor- 
rugation. In order to mix the unperturbed K and K' 
Bloch states, such a profile should break the initial trans- 
lational symmetry in such a way that both points K and 
K' are translated into the origin (F— point) within the 
new "folded" Brillouin zone (BZ). Clearly, this can be 
achieved only if a deformed graphene sheet that (i) still 
has a hexagonal Bravais lattice with the angle 7r/3 be- 
tween two primitive super-periods of equal length, Ai 
and A2, and (ii) its first reciprocal lattice vectors with 
the same length, gi and g2, are specifically commen- 
surate with the vectors K and K' . As to require- 
ment (i), it is easy to list all the hexagonal superlat- 
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tices that are commensurate with the graphene honey- 
comb lattice. Indeed, their periods are Ai = nai -\- ma2, 
X2 = —mai + {n + 771)02, where n and m are arbi- 
trary integers, ai and 02 are the initial graphene pe- 
riods. Accordingly, their reciprocal vectors g defined 
through the corresponding graphene vectors G are: — 
aGi+l3G2, 92 = 1G1+SG2, where a = (n+m)/ f{n, m), 
/3 = m/ f{n,ni), 7 = —m/f(n,m), S = n/ f{n,m); and 
/(ri, m) — + nm + m^. 

From all possible hexagonal superlattices discussed 
above, not all satisfy the condition (ii), but only those 
where 2n -|- m is the multiple of three [3^. Interest- 
ingly, in ID case of (n, m) carbon nanotubes this is a 
condition of their gapless (metallic) behavior, while in 
the present 2D cases it facilitates opening up the gap, 
as we discuss below. Such structures, capable of map- 
ping the graphene K and K' points to the F— point, can 
be divided in three classes: zigzag, armchair, and chiral, 
depending on the directions of A. For the zigzag super- 
lattices, n = 3N, m = 0, so that the periods Ai = 3Nai 
and A2 — 3Na2 are parallel to two different zigzag direc- 
tions. This type of superstructures can also be referred 
to as 3A'' X 3iV, where N is an integer. In the case of 
armchair-edged corrugations, n — m — N. Here, the 
super-periods are rotated by 30° relative to the graphene 
vectors and 02. Since the length of the super-periods 
is larger than that in graphene by a factor of \/3N, the 
structures themselves can be labeled as V3N x V3N. 
And, finally, the chiral superlattices can also be referred 
to as ^ f{n,m) x ^ f{n,m). 

It is worth noticing that the -s/S x superstruc- 
ture breaks the initial translation symmetry in the same 
way as so-called Kekule distortion does, where the in- 
plane atomic displacements can be viewed as a frozen 
symmetrized combination of transverse optical phonon 
modes at i^-points [201 |3T] . While the Kekule distortion 
does not break the point group symmetry of the system, 
Dg/i, this is not the case for all out- of plane 2D corruga- 
tions considered here. Indeed, introducing sine-like cor- 
rugations kills at least the horizontal mirror plane, c/j, 
and the resulting point group is inevitably lower than 
Dg/i. (actually it can be one of its subgroups: Cg^, Dg, 
Ce, Sg, D3d, C3., D3 or C3 [32]). 

All three classes of periodic ripples, in their simplest 
form, can be modeled by the out-of-plane atomic dis- 
placements Uz{r) 



A 



(1) 



where r are the undistortcd lattice positions, A is the 
amplitude, g^, g2, g^ — —gi — ^2 ^^'^ the first supercell 
equal-length reciprocal lattice vectors, each rotated by 
27r/3 with respect to one another, rg is the origin of the 
coordinate system, and (f some phase. Below, we will 
characterize the amplitude A via a dimensionless ratio 



TABLE I: Point groups of 2D corrugated structures as a func- 
tion of ro and (p chosen in Eq. ([T]). Note that the displace- 
ments liz 's ([1]) are identically zero for the = ^/3 x ^/3 struc- 
ture when ro is in the center of a hexagon and (p = 0, it/2, so 
that the structures of this family appear only for N > 1. 
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A/\, where A = |Ai| = IA2I = a^J f{n,m) is the length 
of a super-period; a = |ai| = \a2\- Note that the sum of 
two exponents in Eq. ([T]) reduces to the cosine function 
when iy9 = 0, and to the sine function when (p = 7r/2. 

The displacement field (r), Eq. ([T]), leads to a strain 
field (quadratic in Uz, since the linear terms vanish) 



1 duz 

2 dxi 



duz 
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which is also periodic and has non-zero Fourier compo- 
nents of the type g^ -\- g^ . Since such a strain modulates 
the interatomic hopping integrals, the electrons of the su- 
perlattice should "feel" strongly not only the components 
of the potential with the smallest vectors g^, but also 
those with ztg^ ± g^. This, in turn, should promote the 
gap opening due to enhancing of the second and higher 
order diffraction effects. 

The resulting point symmetry of a superstructure im- 
posed by Eq. ([I]) depends on a choice of the origin 
and phase tp. Table [T] lists, for example, the symmetries 
corresponding to two choices of Tq - cither in the center 
of a hexagon or on a carbon atom - and two choices of 
ip: and 7r/2. Although all the superstructures listed 
in Table |l] have identical in shape (hexagonal) BZs, the 
irreducible part of the latter can differ in going from one 
structure to another. Thus, the irreducible part for the 
Cgi, and D^d superlattices is geometrically identical to 
that for graphene, whereas for Cg, Sg, C^y and D3 is 
twice and for C3 four time larger than that for graphene 
[32] . For the Cgi, and T>3d structures, we denote the high- 
symmetry points at the corners of the superlattice BZ as 
Kg, and at centers of the edges as M^, Fig. [T] 

To illustrate the suggested concept, we have performed 
band structure calculations for the corrugated graphene 
sheets using density functional theory implemented in 
the ABINIT package [33] (for computational details, see 
Supplemental Material ??.) Among the groups listed in 
Table [l| only Cgt,, D3d, Cg and Sg contain the symmetry 
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elements that interchange the A and B sublattices. For 
such groups, the gap opens up solely due to inter-valley 
mixing: the distortion breaks the 4-fold degeneracy at F 
point of a flat graphene into two 2-fold degeneracies. In 
the case of armchair Get, and 3x3 superlattices, 

the gap reaches substantial values of 0.33 and 0.30 eV, 
respectively, as the A/X increases from to 0.06 (Fig. [T]). 
With the constraint A/X — 0.06, these values are reduced 
by a factor of three in passing to the armchair 2^/3 x 
2\/3 superlattices having the same symmetries. They are 
further reduced by a factor 2 — 3 in going to the chiral 
Cg and Se x ^/21 structures corresponding to n = 1 
and m = 4. Clearly, this progressive narrowing of the 
gap is due to increasing period A, because the ratio A/X 
remains the same. 
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FIG. 1: (color online) Gap opening in the zigzag Cei, and 
Dad 3x3 superlattices (F — Ks direction) as the ratio A/X is 
increased from to 0.06. The upper panel shows the periods 
of the superlattices and their irreducible part of the BZ. 

The corrugations associated with the symmetries C^^ , 
D3 and C3 (Table |l]) break not only the translational but 
also the sublattice symmetry. In such an event, both 
the inter-valley and intra-valley mixing effects contribute 
to the reconstruction of the band structure in the vicin- 
ity of the Fermi level. In the case of armchair C^^ and 
D3 1/3 X -s/S structures, for example, the four-fold de- 
generate state at F splits into two single levels and a 
two-fold degenerate state. Fig. [2j The 2-fold degenerate 
state is above the Fermi level, so that the gap opens up 
and becomes indirect when the A/X exceeds some criti- 
cal value (as in Fig. [2]). Although in the zigzag C^y 3x3 
superlattices the 4-fold degenerate state also splits into 




r-g/2 

FIG. 2: (color online) Same as in Fig. 1, but for the armchair 
Ca^ and Da VS x -\/3 superlattices in the F — gi/2 direction. 



three states, here the two-fold degenerate state lies be- 
tween the two single levels and the gap does not open up 
at all, Fig. |3] This result is quite unexpected, it shows 
that breaking the inversion symmetry is not a sufficient 
condition for the gap opening. Note that the principal 
possibility for graphene lattice to have zero energy gap in 
the absence of inversion symmetry was discussed earlier 
in Ref. [33]. 

It is interesting that in going over from x \/3 to 
(3\/3) X (3\/3) structures with the same point symme- 
tries {Csv and D3), the character of splitting of the 4-fold 
degenerate state remains but only for very small corru- 
gation parameters A/X < 0.005. As the A/X is further 
increased, the 2-fold degeneracy is lifted and the struc- 
tures exhibits four single levels. We should stress that 
if only the lowest order effects of the corrugations were 
considered (quadratic in A/X), the splitting of the four- 
degenerate level would follow the ratios obtained from 
symmetry considerations. In the present calculations, 
however, the effects of corrugation are included up to 
an infinite order, so that the character of splitting is not 
in one-to-one correspondence with the symmetry of de- 
formed structures. The 4-fold degeneracy is completely 
lifted at any A/ A for the chiral superlattices with C3 sym- 
metries, as exemplified by the C3 V2I x \/21 structure 
in Fig. |4] Here, a modest gap on the order of 0.1 eV 
opens up as A/X reaches 0.06. Note that once the gap 
is opened at F, it does not necessarily remains the direct 
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FIG. 3: (color online) Same as in Fig.l, but for the C31, su- 
perlattices in the F — {g-^ + g^)!'^ direction. 
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FIG. 4: (color online) Same as in Fig.l, but for the chiral C3 
\/2T X \/2T superlattices in the F — gi/2 direction. 



gap, the conduction band minimum may move around in 
the BZ, and the gap may become indirect, as it is seen 
from the Fig. |4]for the case ip — n/2. The above results 
provide exhaustive classification of possible scenarios for 
gap opening in graphene subject to the two-dimensional 
corrugation. It would be interesting to check the above 
predictions in experiments with various substrates pro- 
ducing planar corrugations of a top graphene layer. 

Supplement: Calculation Method. The DPT cal- 
culations for both flat and corrugated graphene have been 
performed using ABINIT [M] within the local density 
(LDA) approximation. Approximately the same k-point 
density for Monkhorst-Pack fc-point grid [3S] was used for 
all superstructures; roughly it corresponds to a 16x16x1 
grid for a flat graphene with 2 atoms per unit cell. The 
sheets have been simulated by a slab-supercell approach 
with the inter-planar distances of SOas to ensure negli- 
gible wave function overlap between the replica sheets. 
The plane waves pseudopotentials have been chosen in 
the form of Troullier-Martins [36] where carbon 2s and 
2p electrons have been considered as valence states. Por 
the plane-wave expansion of the valence and conduction 
band wave- functions, a cutoff energy was chosen to be 80 
Ry. 
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